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the null results of the other experiments. This type of dark matter can be captured 
by the Earth and shield detectors because it is self-interacting. This effect will lead 
to a diurnal modulation in dark matter detectors. We estimate the size of this effect 
for dark matter detectors in various locations. For a detector located in the northern 
hemisphere, this effect is expected to peak in April and can be detected for opti- 
mistic parameter choices. The diurnal variation is expected to be much larger for 
detectors located in the southern hemisphere. In particular, if the CoGeNT detector 
were moved to e.g. Sierra Grande, Argentina then a 5a dark matter discovery would 
be possible in around 30 days of operation. 
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1 Introduction 



Mirror [1] and more generic hidden sector dark matter models provide a consistent 
explanation [2, 3, 4] of the DAMA annual modulation signal[5], as well as the low 
energy energy excess seen in the CoGeNT[6] and CRESST-II[7] experiments. This 
explanation is consistent with the constraints from higher threshold experiments, 
such as CDMS[8] and XENON100[9] when reasonable systematic uncertainties in 
energy scale are considered 2 . 

A distinctive feature of mirror and more generic hidden sector dark matter models 
is that the dark matter particles are self-interacting. This means that dark matter 
captured by the Earth can effectively block the halo dark matter wind. This effect 
will also lead to a diurnal modulation of detector dark matter signals. The purpose 
of this paper is to study this effect. Although we focus particular attention on 
mirror and more generic hidden sector dark matter models with unbroken U(l)' 
interaction, our study would be of relevance to many other self-interacting dark 
matter models[13]. 

Diurnal modulation of a dark matter signal can also occur due to elastic scatter- 
ing of dark matter on the constituent nuclei of the Earth[14]. In that case, the effect 
is typically small unless the dark matter abundance is much less than expected. We 
will show that the diurnal effect arising from the self- interactions of mirror/hidden 
sector particles can be much larger. In particular, for a detector located in the 
southern hemisphere, the diurnal modulation effect is anticipated to be maximal, 
with a rate suppression varying between zero and near 100% for a detector located 
in e.g. Sierra Grande, Argentina. In fact, a dark matter experiment was conducted 
there in the 1990's[15]. Unfortunately, the energy threshold was too high to detect a 
significant rate of dark matter interactions in that experiment. However, CoGeNT's 
low threshold detector is highly portable and could easily be moved from the Soudan 
laboratory to e.g. Sierra Grande. This would allow the diurnal modulation to be 
detected with 5a C.L. in around 30 days of operation. 

2 Mirror & hidden sector dark matter 

Mirrorfl] and the more generic hidden sector dark matter models[2, 16] assume 
the existence of a hidden sector which contains an unbroken U(l)' gauge interac- 
tion which is mixed with the standard U(l)y via renormalizable kinetic mixing 

2 There are also lower threshold analysis by the XENON10[10] and CDMS collaborations [11]. 
However it has been argued[12] that neither analysis can exclude light dark matter when systematic 
uncertainties are properly taken into account. 
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interaction: [17, 1] 



r . = f^F 1 (1) 



where is the standard U(l)y gauge boson field strength tensor, and F' is the 
field strength tensor for the hidden sector U(l)'. This interaction enables hidden 
sector U(l)' charged particles (of charge Qe) to couple to ordinary photons with 
electric charge Qe'e = ee[18]. [In the case of mirror dark matter, e' = ej cos6^]. 
We consider the case where the hidden sector contains two (or more) stable U(l)' 
charged dark matter particles, F± and F 2 with masses m 1 and m 2 . Under the 
standard assumptions of a dark halo forming an isothermal sphere, the condition of 
hydrostatic equilibrium relates the halo temperature of the particles to the galactic 
rotational velocity, v rot ~ 240 km/s: 

T = \mv 2 rot (2) 

where m = nF i mi + nF 2 m2 j g ^ e mean mass f the particles in the galactic halo. We 

have assumed that the self interactions mediated by the unbroken U(l)' gauge in- 
teractions are sufficiently strong so that they thermalize the hidden sector particles, 
F\ and F 2 . The interaction length is typically much less than a parsec[19] and 
the dark matter particles form a pressure-supported halo. The dark matter parti- 
cles are then described by a Maxwellian distribution with fi(v) = exp(—E/T) = 
exp(—\miV 2 /T) = exp[— v 2 /vq{i)\ where 




(3) 

With the assumptions that m 2 » mi and that the abundance of F 2 is much less 
than Fi, we have that Vq(F 2 ) v 2 0t . The narrow velocity dispersion (recall a 2 = 
3fg/2) can greatly reduce the rate of dark matter interactions in higher threshold 
experiments such as CDMS[8] and XENON100[9] whilst still explaining the signals 
in the lower threshold DAMA and CoGeNT experiments. 

Mirror dark matter refers to the specific case where the hidden sector is an 
exact copy of the standard model sectorfl] (for a review and more complete list of 
references see ref.[20]) 3 . In that case a spectrum of dark matter particles of known 
masses are predicted: e', H', He', O', Fe',.- (with m e i = m e ,mH> = rn H , etc). 
The galactic halo is then presumed to be composed predominately of a spherically 



3 Notc that successful big bang nucleosynthesis and successful large scale structure requires 
effectively asymmetric initial conditions in the early Universe, T" <C T and ny /m ~ 5. See ref.[21] 
for further discussions. 
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distributed self interacting mirror particle plasma comprising these particles [22]. 
Ordinary and mirror particles interact with each other via kinetic mixing of the 
U(l)y and its mirror counterpart. 

Both mirror and the more generic hidden sector dark matter models can explain 
the direct detection experiments [2, 3, 4, 23, 24]. These scenarios involve kinetic mix- 
ing induced elastic (Rutherford) scattering of the dark matter particles off target 
nuclei. In the mirror dark matter case, the DAMA/CoGeNT signals are assumed 
to arise from the scattering of the dominant mirror metal component, A', off tar- 
get nuclei. [The He' and H 7 components are too light to give a signal above the 
DAMA/CoGeNT energy threshold]. In the case of generic hidden sector models the 
role of A' is played by F 2 . Such elastic scattering can explain the normalization 
and energy dependence of the DAMA annual modulation amplitude and also the 
CoGeNT and CRESST-II spectrum consistently with the null results of the other 
experiments, and yields a measurement of ey/^ an d mA'/mF 2 - In the case of mirror 
dark matter the favoured parameter space is [4, 3]: 

e\/& ~ few x 10_1 °> 
mA ' -s 16-56 (4) 



where £4' = tia'Ttia' /(0.3 GeV/cm 3 ) is the halo mass fraction of the species A' and 
m p is the proton mass. In the case of the generic hidden sector model, e is about an 
order of magnitude larger. [The difference arises because in our notation the kinetic 
mixing induces an electric charge eZ'e for mirror dark matter and ee for the generic 
hidden sector case]. 

As briefly discussed in Eq.(3) above, the velocity dispersion of the halo A'/F 2 
particles depends on the parameter fh. In the mirror model with isothermal halo it 
can be estimated[25] to be around 1.1 GeV although it might vary somewhat from 
this value. For example, if the e' have a lower temperature than the mirror nuclei 
in the halo then fh can be somewhat larger, with an upper limit of around 3 GeV. 
Although realistically it is unlikely that fh could be this large in the mirror model, 
we consider this value in our numerical work as an extremely optimistic case which 
might also be relevant to more generic hidden sector models where the parameter 
fh is less constrained. 



3 Distribution of mirror dark matter captured by 
the Earth and its shielding radius - Rq 

In the following sections we restrict our attention to the case of mirror dark matter 
for definiteness. The more generic hidden sector scenario is completely analogous. 
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Mirror particles will occasionally be captured by the Earth and accumulate in 
the Earth's core. Once a significant population of mirror particles have been accu- 
mulated, mirror dark matter will be captured at a rate: 

^ - «I%fA> (5) 

where Ja' ~ v ro t £' A i 0,3 ( ^ V / cm is the flux of A' mirror particles hitting the Earth. 
Here Rq is the maximum distance from the Earth's center for which dark matter can 
be captured due to self interactions. We will show in a moment that this distance 
is (currently) of order 4, 000 km for the example where m A i 22m p . Thus we 
anticipate that around 



N ~ J 7iR 2 f A ,dt 



~ 11)39 (wi) W 

A' particles will be captured during the five billion year history of the Earth. This 
is many orders of magnitude within the geophysical limits [26]. 

A mirror particle, once trapped in the Earth, will lose energy rapidly due to in- 
teractions with the ordinary matter. The Rutherford cross section increases rapidly 
as the relative velocity decreases, da/dfl oc 1/f 4 , and thus one expects the captured 
mirror particles to rapidly thermalize with the ordinary matter. Their distribution 
can then be obtained from the condition of hydrostatic equilibrium: 



dp 




dr 


-PA' 9 


dn A i 


n A > 


dr 


T 



dT 

rriA'g + ~r 
dr 



(7) 



where we have used P = n A iT and p& = m A in A i. Note that we have implicitly 
assumed that captured mirror nuclei and mirror electrons combine into mirror atoms. 
The gravitational acceleration, g, is given in terms of the mass density p, 

g = % [ r p47tr> 2 dr' . (8) 
r z Jo 

Note that both g and T depend on r. We model the Earth's density and temperature 
as given in figures la, lb. The boundaries in figure la correspond to the inner core, 
core, and mantle. In each of these regions we use a linear approximation for the 
density adapted from the Preliminary Reference Earth Model [27]. 
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Figure la: Earth's mass density [g/cm 3 ] versus radial distance. 
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Figure lb: Earth's temperature profile [°K] versus radial distance. 
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Eq.(7) can be numerically solved for nA'(f), taking into account the dependence 
of g and T on r. The result is shown in figure 2. If the total number of accumulated 
A' particles is around 10 39 , as indicated in Eq.(6), then we find that the central 
number density is of order n^'(0) ~ 10 14 cm~ 3 . 




1000 2000 3000 4000 5000 6000 

Radial distance (km) 

Figure 2: The distribution within the Earth, riA'(r)/nA'(0), of captured mirror particles 
of mass niA' = 22m p (solid line), tua' = 30m p (dashed line), mA> = 14m p (dotted line) 
and niA' = ^m p (dashed-dotted line). 



Figure 2 suggests that any captured He' particles will quickly escape the Earth. 
Indeed, the He' distribution is quite spread out, and particles in the tail of their 
Maxwellian velocity distribution can have velocities much greater than the Earth's 
escape velocity. One can easily show[23] that e.g. He' particles with velocities ~ 20 
km/s can travel ~ 100 km through ordinary matter without significant energy loss 
and escape the Earth's gravity (assuming here e ~ 10~ 9 ). Heavier mirror particles on 
the other hand are much more tightly confined within the Earth and have negligible 
probability of escaping Earth's gravity. 
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The captured mirror dark matter in the Earth can potentially shield a detector 
from halo dark matter particles if they originate from a direction which takes them 
through the Earth. Whether or not a halo dark matter particle will get shielded 
will depend on how close its trajectory takes it to the center of the Earth, where 
the density of Earth bound dark matter particles is greatest. Let us define r$ as the 
distance of closest approach to the Earth's center of a given particle trajectory. Let 
us now work out the maximum value of tq, which we call Rq, for which halo dark 
matter will be captured due to self interactions with Earth bound dark matter. 

When a dark matter particle, A', from the halo passes through the dark matter 
accumulated in the Earth it will lose energy. The rate of energy loss is given by: 

d E ' , x f E % ax da „ 7rn 

where x is the path of the dark matter particle. If the density of Earth bound dark 
matter particles is great enough, the dark matter particle A' will lose nearly all of 
its energy from self-interactions and become trapped within the Earth. The relevant 
cross-section for Rutherford scattering of a mirror nuclei with incoming velocity v 
on the mirror nuclei trapped within the Earth (here both assumed to have atomic 
number Z') and neglecting the form factors 4 is given by 

da A (10) 



dE R E 2 R v 2 
where 

A S *^. (11) 

With this cross-section, Eq.(9) indicates that dark matter particles will become 
trapped within the Earth if their trajectories are such that: 



J nA'(x)dx 



E in 



m A i\ log 



Z (|J)Vcm- 2 , (12) 



where E' in ~ \ m A'vl ot is the initial energy of the dark matter particle. The quantity 
£max _ ig maximum kinematically allowed E R , while E R im ~ l/(2rfm,A') (n 
is the Bohr radius) is set by the scale at which atomic screening becomes important. 



4 We employ natural units where h = c = 1. 
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Approximating the trajectories of dark matter particles by straight lines, and 
defining the co-ordinate q as the distance along the trajectory with r equal to the 
distance of closest approach to the Earth's center, we have: 



Eq.(12,13) can be numerically solved for r assuming a given value. We find 
numerically that dark matter particles will get captured by self interactions provided 
their trajectories satisfy ro ~ Ro, with Rq ps 4,000 ± 700 km, for m.A' /m p = 22 ± 
8 5 . Uncertainties in the p,T profile of the Earth give an (additional) estimated 
systematic uncertainty of order 5 — 10% in the determination of Rq for a given 
value. Thus, R could in principle be as large as 5, 500 km. This completes our 
estimation of the 'shielding radius' R . 

4 Diurnal modulation effect due to shielding of 
the halo dark matter wind 

We are now ready to estimate the effect of the shielding of the halo mirror dark 
matter due to its interactions with mirror dark matter trapped within the Earth. 
In principle halo dark matter interactions with ordinary matter can also shield or at 
least reduce the energy of halo dark matter particles[14]. However, for the case of 
mirror dark matter the latter effect is expected to be relatively small [23] especially 
if the kinetic mixing parameter e satisfies the cosmology bound[28] e ~ 10~ 9 . 

The direction of the Earth's motion through the halo, subtends an (average) 
angle ~ 43° with respect to the Earth's spin axis. There is a small annual modulation 
in this angle due to the Earth's rotation around the sun, which can be important 
and will be considered in the following section. Another relevant angle is the angle 
between the direction of the Earth's motion through the halo and the normal vector 
to the Earth's surface at the detector's location. This angle is denoted by ip [Note 
that in ref.[14] this angle was denoted as 6]. The angle ip, depends on the detector 
location and also the time of day, t: 



where ui = 2ir/T d with T d = 1 sidereal day (23.934 hours). In the above equation, the 
+ [— ] sign is relevant for a northern [southern] hemisphere detector, where 9i a titude 
is the north [south] latitude. 

5 Note that even if heavy <~ Fe' component is responsible for the DAMA, CoGeNT and CRESST- 
II signals [4], it is reasonable to expect lighter components do dominate the shielding within the 
Earth. 




(13) 



cos ip = cos Oiautude sin ut sin 43° ± sin 6i atitude cos 43° 



(14) 
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Dark matter particles arriving at a detector originate predominately from a cone 
with axis in the direction of the Earth's motion through the halo and with the 
detector at its apex. The particles within the cone which pass through a distance 
from the Earth's center of r ~ R (R pa 4, 000 km) will be shielded from the 
detector, while particles with r ~ Ro will arrive unhindered by self-interactions. 
Because the direction of the cone of dark matter particles periodically changes due to 
the Earth's rotation, the proportion of dark matter particles shielded also modulates. 
Our task now is to estimate this effect for detectors in various locations. 

The differential interaction rate is given by: 

dE~ R = NtU « } dE R —T- |v|<i " 

= Km* r Mjfflfv (15) 

E Z R J\ v \>v mm (E R ) k\v\ 

where Nt is the number of target atoms per kg of detector, k = (ttvq[A']) 3 ^ 2 is the 
Maxwellian distribution normalization factor and riA> = Pdm£,A> jmA> is the num- 
ber density of the halo mirror nuclei A' at the Earth's location (we take pd m = 
0.3 GeV/cm 3 ). [This should not be confused with the number density of trapped 
mirror nuclei in the Earth, despite the similarity in the notation!]. Here v is the 
velocity of the halo particles relative to the Earth and v# is the velocity of the Earth 
relative to the galactic halo. The halo distribution function in the reference frame 
of the Earth is given by, f A '(v, v E )/k = (irvo[A'})~ 3 ^ 2 exp(— (v + ve) 2 /vq[A']). Note 
that the lower velocity limit, v m i n (En), is given by the kinematic relation: 



(m A + m AI ) 2 ER ^ 



u min \ n 9 ' 

\ 2mAm AI 

The interaction rate depends on the velocity integral in Eq.(15), 

/a'(v,v b ) 3 



/ ee / J IY ; d*v . (17) 

JM>v mln (E R ) fc|v| 

This integral can be modified to incorporate the effect of shielding of the halo dark 
matter by the Earth bound mirror matter, by multiplying the integrand by the 
quantity g(9, 0, ip), where 

g(6,(j),ijj) = if d min < R 

= 1 if d min > Ro ■ (18) 

Here d min is the distance of closest approach to the center of the Earth experienced 
by the particle trajectories. We show in the appendix that d m i n is given by: 

<& in (M,V0 = R 2 E [i-f\BAM^ MW)>o 

= R 2 E ii /(M,VO<0 (19) 
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where 



f(6,(f>,ip) = sin 9 sin ip sin — cos 9 cos ip (20) 

and Re — 6378 km is the radius of the Earth. Thus, the shielding effects can be 
incorporated by replacing I, defined in Eq.(17), with: 



\v\>v min (E R ) 

i r 2 - /■« /a'(v,v b ) 



/[^)] = T J^^- giO^Mv 

JW\>v min (E„) fc|V] 

-#(0,0,-0)|vM|v|#dcos0 . (21) 



1 JO J|v|>D min (S i j) fc 

We are now in a position to evaluate the effect of shielding of mirror dark matter 
expected for experiments at various locations. We compute the quantity 

K = 100 x (l - SipPj (22) 

which gives the percentage rate reduction in the interaction rate due to shielding 
of dark matter in the core of the Earth. In figure 3a we give results for 1Z, with 
Er = 6.7 keV^R for Na target, relevant for the DAMA set up at Gran Sasso with 
north latitude = 42.5°. We consider two reference sets of parameters, the first one 
assumes the 'standard' parameters Rq = 4,000 km and m = 1.1 GeV. The second 
possibility assumes 'optimistic' values R = 5, 500 km and m = 3.0 GeV. Figure 3a 
shows that even with the optimistic parameters, the diurnal signal modulation is 
expected to be relatively small for DAMA/Libra at Gran Sasso. Thus the lack of 
a diurnal modulation in the DAMA/Nal experiment [29] does not significantly con- 
strain any parameter space. Currently DAMA/Libra is running with a lower energy 
threshold which will make them more sensitive to the diurnal effect simply because 
the signal rate is much higher at lower energies, dR/dE R oc 1/E R . Thus, for opti- 
mistic parameters DAMA/Libra might possibly detect a small diurnal modulation. 

The CoGeNT Ge detector has been operating in the Soudan underground lab- 
oratory, which has north latitude of 48°. This is a little further north than that of 
Gran Sasso, so the diurnal modulation effect expected at Soudan is smaller than 
that for DAMA/Libra. However, for the CDEX/Texono Ge detectors[30] located at 
Jin-Ping underground laboratory at north latitude 28°, the effects are expected to 
be larger. In figure 3b we give the results for these detectors at E R = 2 keV^R- 
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Figure 3a: Percentage rate suppression due to the shielding of dark matter in the Earth's 
core versus time, for a Na detector at Gran Sasso. A reference energy of 6.7 keV^R has 
been assumed. Both the captured and interacting dark matter particles are assumed to 
have the mass = 22m p . The solid line assumes 'standard' parameters where the 

shielding radius is taken to be Rq = 4, 000 km and the halo velocity dispersion is given in 
Eq.(3) assuming fh = 1.1 GeV. The dashed line is for the more optimistic case considered 
where Rq = 5, 500 km and fh = 3.0 GeV. 
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Figure 3b: Same as figure 3a except that the Na detector at Gran Sasso is changed to 
CDEX/Texono Ge detector at Jin-Ping underground laboratory (a reference energy of 
2 keVxR is assumed). 



The effects in existing detectors in the northern hemisphere are not very large, 
although there is some possibility of observing the effect with optimistic parameters, 
especially at the Jin-Ping laboratory. Of course, the reason that the shielding effect 
is not so large is simply because all of the existing detectors are located in the 
northern hemisphere. The effect becomes dramatically bigger for a detector located 
in the southern hemisphere because the halo dark matter particles then have a much 
larger probability of passing through the center of the Earth We illustrate this in 
figure 3c which gives results for the CoGeNT detector relocated to the Sierra Grande 
laboratory in Argentina, with south latitude 41.6°. As the figure shows, the variation 
ranges from near zero to 100 %. At Soudan, the CoGeNT spectrum is consistent 
with around 3 dark matter signal events per day per 0.3 kg. We estimate then 
that the diurnal effect at Sierra Grande would be so large that just 30 days of data 
would be sufficient to yield a 5a discovery of dark matter with CoGeNT's ~ 0.3 kg 
detector. 
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Finally, we comment on the proposal to place a Nal detector at the South 
Pole [31]. Interestingly, the South Pole is the only point in the southern hemisphere 
where the diurnal modulation effect is absent. For a detector at the south pole we 
find a constant 32% [85%] rate suppression for the standard [optimistic] parameter 
values. 



5 Annual variation of the diurnal modulation 

The proportion of the Earth which shields dark matter particles has both a diurnal 
variation and an annual variation. The annual variation, which we have hitherto 
neglected, arises because of the variation of the direction of the Earth's motion 
through the halo, due to the Earth's motion around the sun. In this section we 
examine this effect. 
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The direction of the Earth's motion through the halo, subtends an angle 9\(t) 
with respect to the Earth's spin axis. On average this angle is approximately 43°, 
but varies during the year. To evaluate this angle at a given time, t, define a 
right-angled co-ordinate system with the z-axis in the direction of the normal to the 
ecliptic plane, x-axis in the direction of the sun relative to the Earth and y-axis in 
the direction of the Earth's motion around the sun. In this co-ordinate system, for 
an observer at rest with respect to the halo, the Earth's velocity vector is 

v E = i; [-sin3O°cos27r ^~ tl \ - sin 30° sin 2tt ^ ~ fl \ cos 30°] + v @ [0, 1, 0] 

= v Q [- sin30°cos27r ^ - sin 30° sin 2tt ^ ~ ^ +y, cos 30°] (23) 

where t\ = 152.5 days + 0.25 years ~ 244 days, T = 1 year and y = i>®/i>© ~ 
30/240 pa 0.125. The angle, 30° is the angle between the normal of the ecliptic plane 
and the sun's direction of motion through the halo. The Earth's spin axis has the 
direction: 

L spin = [sin 23.5° cos 2tt ^ ~ ^ , sin 23.5° sin 2? ^ ~^ 2 \ cos 23.5°] (24) 

where t 2 ~ 172 days (northern summer solstice). The angle 23.5° is the tilt of the 
Earth's spin axis relative to the normal of the ecliptic plane. The angle 9i(t), is then 
given by 



cos#i = 

~ cos 9~i + y 

where 



cos 6 1 sin 30° sin 2vr-^ — ^ + sin 23.5° sin 2tt^ — ^ 

T T 



(25) 



cos 0i = cos 30° cos 23.5° - sin 30° sin 23.5° cos 2?r ^ 1 (26) 

Thus, we see that the angle 9\{t) is not actually a constant, but has a small variation 
due to the earth's motion around the sun (y ^ 0). 



14 



60 



50 - 



40 - 




03 

S 30 - 
20 - 



10 - 



I i i i i i i i I 

50 100 150 200 250 300 350 

Time (days) 

Figure 4: The angle, 9\{t) versus time. [This is the angle between the Earth's spin axis 
and the Earth's direction of motion through the halo]. 

In figure 4, we plot 9\(i) over its period of one year (where t = corresponds to 
January 1). As the figure shows, we find that 9\ has a maximum of 49.3° (April 25) 
and minimum of 36.3° (October 25). This variation is relatively unimportant for a 
detector in the southern hemisphere 6 . For a detector in the northern hemisphere 
we expect the largest diurnal signal to occur when 9± is a maximum and thus we 
anticipate that the largest diurnal signal will occur in April. The effect of the d\ 
variation can be estimated by repeating the analysis of the previous section, with 43° 
in Eq.(14) replaced with 49.3 and 36.3 [in general 43° is replaced by 0i(t) in Eq.(14)]. 
In figure 5a we plot the percentage rate suppression for the optimistic parameter 
choice for a detector located at Gran Sasso. The figure shows that the variation 
of the diurnal signal during the year is significant. Figure 5b is the corresponding 
results for a detector located at Jin-Ping Underground laboratory. Clearly the March 

6 For the special case of a detector at the South Pole, there is no diurnal signal as noted earlier, 
but there is a rate suppression which annually modulates slightly due to the variation of 9\{t). 
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- May quater provides the best window to detect the diurnal signal for a northern 
hemisphere detector. The CDEX experiment in Jin-Ping, currently running with a 
1 kg Ge PC detector, with plans to upgrade to > 10 kg has an excellent opportunity 
to observe the diurnal signal. The 250 kg DAMA/Libra detector, currently running 
with lower threshold, also has a very good chance to see such a signal in the March- 
May quater. 
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Figure 5a: Percentage rate suppression due to the shielding of dark matter in the Earth's 
core versus time, for a Na detector at Gran Sasso. A reference energy of 6.7 keV^n has 
been assumed. Both the captured and interacting dark matter particles are taken to have 
the reference mass m^ 1 = 22m p . We have assumed the optimistic case with shielding 
radius of Rq = 5, 500 km and the halo velocity dispersion is given in Eq.(3) with fh = 3.0 
GeV. The solid line assumes 6\ = 49.3° (April 25), the dashed line assumes 6\ = 43°, 
(Yearly average) and the dotted line assumes 0\ = 36.3° (October 25). 
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Figure 5b: Same as figure 5a except that the Na detector at Gran Sasso is changed to 
CDEX/Texono Ge detector at Jin-Ping underground laboratory (a reference energy of 
2 keVNR is assumed). 

6 Conclusion 

Mirror and more generic hidden sector dark matter models can simultaneously ex- 
plain the DAMA, CoGeNT and CRESST-II dark matter signals consistently with 
the null results of the other experiments. This type of dark matter can be captured 
by the Earth and shield detectors because it is self-interacting. This effect will lead 
to a diurnal modulation in dark matter detectors. We have estimated the size of 
this effect for dark matter detectors in various locations. 

Current dark matter detectors in the northern hemisphere are not so sensitive 
to the diurnal modulation although there is an exciting possibility that it can be 
detected at Jin-Ping and even in DAMA/Libra for optimistic parameter choices. 
We have also found that the diurnal signal for detectors located in the northern 
hemisphere varies during the year, peaking in April. 
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The situation changes dramatically for detectors located in the southern hemi- 
sphere where the dark matter wind passes directly through the core of the Earth for 
part of the day. In particular, if the CoGeNT detector were moved to e.g. Sierra 
Grande, Argentina then a 5a dark matter discovery would be possible in around 30 
days of operation. 

Appendix: Calculation of d min 

In this appendix we derive Eq.(19). Recall dark matter particles reaching a 
detector originate predominately from a cone centered along the Earth's direction 
of motion through the halo. Recall also that ip is the angle between the normal vector 
to the Earth's surface at the detector's location and the axis of this cone. We assume 
that particles form straight line trajectories. Let us introduce a co-ordinate system 
with the detector at the origin. We use polar co-ordinates with the z — axis aligned 
with the direction of the Earth's motion through the halo. The points along the 
particle trajectory are given by (r sin 9 sin 0, r sin 9 cos 0, r cos 9). In this co-ordinate 
system (suitably rotated in the azimuthal direction) the location of Earth's center is 
sin-0, 0, —R E cos?p), where Re — 6378 km is the Earth's radius. The distance 
between the Earth's center and a point along the particle trajectory is then given 
by 

d 2 = R\ + r 2 + 2R E r(cos 9 cos ip - sin 9 sin ip sin 0) . (27) 

We need to determine d m i n which is the distance of closest approach to the Earth's 
center as we vary r keeping 9, and ip fixed. The value of r where this occurs can 
be determined by minimizing d 2 with respect to r, i.e. from dd 2 /dr = 0. Evaluating 
this value of r and substituting it back into Eq.(27) we obtain Eq.(19) as the solution. 
This completes the outline of the derivation of Eq.(19). 
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